Quantum computation with ultracold atoms in a driven optical lattice 
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We propose a scheme for quantum computation in optical lattices. The qubits are encoded in the 
spacial wavefunction of the atoms such that spin decoherence does not influence the computation. 
Quantum operations are steered by shaking the lattice while qubit addressability can be provided 
with experimentally available techniques of changing the lattice with single-site resolution. Numer- 
ical calculations show possible fidelities above 99% with gate times on the order of milliseconds. 



In the last years tremendous progress has been made in 
controlling and observing ultracold atoms in optical lat- 
tice (OL) potentials [1 . One of the latest developments 
has been the optical detection of atoms with single-site 
resolution in lattices of increasingly smaller periodicity 
[2H5]. This technological advancement allowed, e.g., for 
a direct observation of the superfluid to Mott-insulator 
transition [6J . Along with these detection schemes comes 
the possibility to control the lattice potential with single- 
site resolution. Lately, this has been used to manipulate 
the spin of single atoms in an OL [7]. 

One of the possibly most important applications of 
these techniques is the implementation of a quantum 
computer that would dramatically improve the compu- 
tational power for particular tasks. Compared to other 
possible candidate systems for the implementation of a 
quantum computer, neutral atoms in OLs have the ad- 
vantage of a natural scalability to a large number of 
atoms encoding qubits and a weak coupling to the envi- 
ronment leading to long decoherence times ^ . Single-site 
addressability of large qubit systems may be one of the 
last milestones on the way to an OL quantum computer. 

Many schemes proposed for quantum computation in 
OLs rely on encoding the qubits by atomic spin states j9r 
[T^ . Although single-qubit operations and collective 
two-qubit operations '13] have been demonstrated, spin 
qubit states are generally disturbed by external mag- 
netic fields that lead to their decoherence. This source 
of decoherence can be avoided by encoding the qubits 
in the spacial wave function of the ground and first ex- 
cited state of Bosonic atoms localized at single sites of 
an OL [121 [HI [15]. In this Letter we show that a pe- 
riodic modulation of the lattice position, i.e. a shaking 
of the lattice suffices to perform all quantum operations 
needed for quantum computation. The shaking of the 
OL can be studied in the context of Floquet theory show- 
ing that it effectively changes the hopping parameter of 
the system [TB] . This effect has been verified experimen- 
tally, while revealing that the shaking also drives transi- 
tions from the first to the second Bloch band [T7] . Qubit 
addressability is provided by manipulating the OL with 
single-site resolution. This enables the selective change 
of the energy spacing of specific sites and the driving 
of local transitions. Additional z rotations induced by 



the lattice manipulation can be cancelled by refocusing 
schemes known from NMR quantum control [TS] . 

In the following we consider a one-dimensional (Id) 
OL in X direction with a potential Vl = Vosin^(fca;) as- 
suming tight confinement in the transversal directions. 
Lengths are given in units of 1/fc and energies in units 
of the recoil energy Er = h'^k'^/{2m). In this system of 
units the lattice spacing is given as d = tt and the oscil- 
lator energy of the harmonic approximation of a single 
lattice site as huL — 2-\/K)- Ultracold Bosons in OLs are 
regularly described by the Bose-Hubbard model which 
uses the Wannier basis of the first Bloch band to for- 
mulate the Hamiltonian in second quantization [19) . In 
order to consider excitations to higher Bloch bands the 
basis has to be extended to include Wannier functions 
Wi^b{x) of sites z = 1 ... TV and several bands b = 1 . . .n. 
Generally, atoms in different Bloch bands are coupled by 
the interaction Vint(a;i, 0:2) = gS{xi — X2)- By means of 
a Feshbach resonance or strong transversal confinement 
one can make g sufficiently large to form a Mott insula- 
tor state in the first two Bloch bands with one atom per 
site. A lattice with unit filling then realizes a quantum 
register with |0) encoded by an atom in the excited state, 
and |1) encoded by one in the ground state [see Fig. [l] 
a)]. In the following 
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shall denote a Fock state in the extended Wannier basis 
with occupation number n^ j of band i and site j. If the 
coupling between the bands and sites is small and no 
avoided crossing appears all important eigenstates can 
be characterized by their dominantly contributing Fock 
state. In this notation the qubits state |001) is encoded 
by the eigenstate with maximal overlap to ||ooi)- 

In the co-moving frame of the lattice a shaking is de- 
scribed by Vsh = fshX cos{uJsht — where Wsh is the 
frequency, if the phase, and /sh(i) the force amplitude 
of shaking [50]. This perturbation leads to a coupling 
between Bloch bands with symmetric and antisymmetric 
Wannier functions [see Fig.[l]b)]. We note in passing that 
transitions between the ground and excited states can 
also be achieved by Raman transitions [21]. By means 
of Vsh one is generally able to drive transitions between 
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FIG. 1. (Color online) a) Qubit register with atoms (circles) 
in the ground (||i') = |1)) or excited state (j|o) = |0)) of each 
lattice site, b) Shaking of the lattice couples predominantly 
ground and excited states in each lattice site (indicated by 
arrows), c) By shining a laser with waist on the order of the 
lattice spacing (red shading) onto one site a local transition 
between |0) and |1) can be driven by shaking the lattice. 

different eigenstates of the lattice Hamiltonian Hq. For 
this consider a Hamiltonian ?^(i) = Ho + cos{ujt — ip), 
where Hq has two eigenstates |^i), \ip2) with eigenener- 
gies hioi and hLj2 , respectively. For u = 102 — uji the full 
time evolution operator in the rotating frame is given in 
terms of the Pauli matrices [a^^Oy] as 

lJ{t) = exp [cos((^)(7a; - sm{ip)ay]^ (2) 

where fi^ = A|(i/'i|V^|V'2)|/^ is the Rabi frequency [22]. 
This enables arbitrary rotations about any axis in the xy 
plane of the system's Bloch sphere. A z rotation may be 
driven by combining rotations about the x and y axes. 

Driving a single-qubit operation at a certain lattice 
site necessitates that the energy difference between |0) 
and |1) differs from that of other lattice sites. This can 
be achieved by shining a laser with a waist on the or- 
der of the lattice spacing perpendicular to the lattice, as 
was accomplished experimentally in [7]. The perturba- 
tion by this laser can be assumed to have the Gaussian 
form Vcauss = -7exp[(a: - a;o)^/(2(T^)]. Depending on 
whether the atoms are strong or weak-field seekers one 
has 7 > or 7 < and the additional intensity will en- 
large or reduce the energy difference h(uj2 — wi) between 
ground and excited state [see Fig.[l]c)]. By shaking with 
a frequency Wsh — LO2 — UJ1 one is then able to drive single- 
qubit rotations on the marked lattice site. 

The detuning of the energy levels leads to an addi- 
tional rotation about the z axis in the Bloch sphere of 
the marked and also of neighboring qubits. Moreover, off- 
resonant shaking of a lattice site can also induce a slight 
phase shift [H] leading to a general dephasing of the 
qubits. These effects lead predominantly to additional 
terms in the Hamiltonian of the form W = J2j o-j ^z'\ 



where the Pauli z operator ctz acts on qubit j. If the a^'s 
are known one can account for the additional z rotations 
within the quantum calculation. However, one can also 
cancel these rotations by applying a scheme similar to the 
refocussing technique well known from NMR quantum 
control [18]. Suppose an x rotation xj^^ = exp(j(/)(7i*''/2) 
is to be driven on a qubit at site s. The time 2t of the op- 
eration is chosen such that exp(— irasfri'V^) = 1- Then, 
inserting a global tt rotation X|' = exp(i7r/2 J^j o'i"'^) be- 
fore and in the middle of the operation does not perturb 

^ (s) 

the ' rotation but all z rotations are cancelled since 

X|' has the property that xfe-'^^^'^X^^e-'*^/'' = 1. 

Together with single-qubit operations, the ability to 
drive a controlled rotation (GROT) between adjacent lat- 
tice sites completes a universal gate set [23]. A GROT 
rotates one qubit (the target qubit) if and only if another 
qubit (the control qubit) is in state The strategy 
to perform this operation between neighboring qubits is 
to deform the lattice such that a repulsively bound state 
lloo) [21] comes into resonance with the state |00) = ||oo). 
The coupling between the lattice sites leads to an avoided 
crossing in the energy spectrum (see Fig.[2|. If we iden- 
tify the left qubit with the control qubit, any rotation 
like that of Eq. ([2| on the right target qubit becomes off- 
resonant and is inhibited iff the control qubit is initially 
in the excited state (|0)). 
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FIG. 2. (Color online) A laser (red shading in right image) is 
positioned slightly offset from the middle between two lattice 
sites so that both lattice sites are coupled and the energy of 
the right lattice site is lower than the one of the left site. At 
a certain laser intensity the energy of the repulsively bound 
state ||oo) and the state \\oo) form an avoided crossing (left 
image). It is now possible to shake the lattice with frequency 
Q resonantly to the transition |10> ^ |11> {\\%) ||??)) 

while the transition |00) |01) (||oo) ^ ||oi)) is off-resonant. 
This enables rotations on the right target qubit conditioned 
by the left control qubit. 

In order to validate the proposal for single and two- 
qubit operations by a numerical study a lattice with unit 
filling described by a Wannier basis of the first three 
Bloch bands with periodic boundary conditions is con- 
sidered. It is assumed that the lattice is sufficiently deep 
so that, as usual within the Bose-Hubbard model, only 
next-neighbor hopping and onsite interaction need to be 
considered. For the single-qubit rotation, a system of 
two lattice sites suffices to estimate the influence of the 
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operation on the remaining system and vice versa. For 
studying the two-qubit operation a third site has to be 
added. The third Bloch band is included to study pos- 
sible excitations of atoms out of the qubit basis. In the 
following a lattice depth Vq = 2.7hujL — 29.16 is consid- 
ered. For this lattice an interaction strength g = 1.87 
suffices to form a Mott insulator so that all qubit states 
are eigenstates of the Hamiltonian. 

As an example for a single-qubit operation the NOT 
operation on the right site of a two-well lattice is con- 
sidered. As stated above the dephasing caused by terms 
of the form W — ai&i^^ + a^c^f^ is inhibited by two 
global X|' rotations. However, the narrow-waist laser 
can also lead to a weak coupling fcai^' • (j'z^ which can 
only be cancelled by more complex sequences of refo- 
cussing pulses. To suppress this term during the op- 
eration the lattice depth Vq is temporarily enlarged to 
Vo -|- SVq before shining in the narrow-waist laser. Both 
perturbations are switched on and off adiabatically and 
are sufficiently small so that couplings to bands above 
the third Bloch band are negligible. While the pertur- 
bations are active V^sh drives resonantly in the adiabatic 
basis the rotation IJ — exp(iSl7?i(7a;/2) on the right site, 
which for — -k is the NOT operation up to a global 
phase. Fig. [3] shows numerical simulations for the NOT 



pulse 1 pulse 2 pulse 1 pulse 2 




SO 100 150 200 250 300 

t \fijEA 



FIG. 3. (Color online) Admixtures of different eigenstates 
for the initial states and jOl) during a NOT operation 
on the right qubit of a two-well lattice. Pulse 1 drives an 
rotation. Pulse 2 starts by linearly increasing the lattice depth 
by (SVo = 0.17Vo during a ramping time of 5.9ft/-Er. Then the 
Gaussian beam with waist a = ■k/2 centered at the site of 
the target qubit is linearly ramped during a time of OAh/E,. 
to the strength 7 = 2.62. After a 7r/2 rotation (VNOT) on 
the right qubit driven by shaking the lattice resonantly to the 
energy difference between jll) and jlO), SVo and 7 are ramped 
off in reversed order. The simulations show a fidelity of 99.7% 
for a gate time of 300h/Er. 

operation including two A"f' operations for refocussing. 



For the chosen parameter set (see caption of Fig. |3]) the 
influence of the qubit state on the left side is negligible 
and an average fidelity [25] of 99.7% is reached. Without 
refocussing (not shown) the fidelity would be only 33.2% 
due to dephasing. 

For the GROT operation a system of three lattice sites 
is considered. The left site acts as the control qubit of an 

' (2) 

Xir rotation on the central target qubit. Up to a local 
phase this operation is equivalent to a controUed-NOT 
operation. Due to the periodic boundary conditions the 
right site is coupled to the control and the target qubit. 
Example results and parameters of the numerical simula- 
tions of the GROT are shown in Fig. |4] For a gate time of 
400h/Er an average fidelity of 99.4% is reached. For the 
evaluation of the fidelity the individual z rotations and 
those z rotations coupled to the control qubit are ne- 
glected. This is possible since equipped with the global 
and local NOT operations these z rotation can be 
easily cancelled by using a refocusing scheme. Opposing 
to the NOT operation an important source of infidelity 
of the GROT stems from the leakage to states out of 
the qubit basis. Averaged over the computational ba- 
sis the leakage probability after the operation is 0.27%. 
To diminish this one can either use techniques of leak- 
age elimination or one has to choose a deeper lattice. 
This leads to weaker coupling between the lattice sites 
and thus to longer gate times. Increasing, e.g., the lat- 
tice depth Vq from 2.7hujL to 2.8fiujL and the gate time 
to 54:0h/Er reduces the leakage probability to 0.15% and 
increases the fidelity by another 0.06%. 

Independently of the chosen lattice depth, the time 
scale h/Er is a crucial system parameter for the speed 
and accuracy of the manipulations. In general a large 
recoil energy Er — allows shorter time scales. 

For example, the NOT operation for a *^Rb system in a 
d = 500 nm OL (Er = 1.5 x lO'^o J) lasts 21 ms while the 
same operation would take 0.6 ms for ^Li in a d = 300 nm 
lattice {Er = 5.2 x lO'^^ J). 

In order for the operations to be robust any relative 
energy shift of the manipulated qubit states by some 
external perturbation must be small compared to the 
energy scale of the Rabi oscillations fi-Jl/j which in our 
case is about 0.1 E'^. Since the energy differences are 
mainly influenced by the lattice potential itself, its un- 
controlled perturbation can severely reduce the fidelity of 
the operation. Gonsidering the GROT operation, which 
is more sensitive to perturbations than single-qubit oper- 
ations, only a lattice laser intensity which is controlled on 
the 10^'* level leads to negligible fidelity reduction (see 
Tab. [l]). On the other hand, perturbations of other pa- 
rameters such as the intensity of the gaussian laser beam 
or the interaction strength are much less severe. 

The last important ingredient for quantum computa- 
tion is the possibility to read out qubit states with high 
fidelity. This can be done by removing atoms in the ex- 
cited state |0) from the lattice and determining subse- 
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TABLE I. Average fidelity of the CROT operation presented 
in Fig. |4] for uncontrolled errors of the interaction strength g, 
the strength of the gaussian laser 7, and the lattice depth Vq. 



Parameter 


g 


7 


Vo 


Error 


0.1% 


1% 


0.1% 


1% 


0.01% 


0.1% 


1% 


Fidelity 


99.4% 


98.1% 


98.9% 


97.3% 


99.2% 


90.9% 


31.4% 
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ing a quantum gas '53^ supporting the practicability of 
vibrationally encoded qubits and the proposed read-out 
scheme. 




FIG. 5. Energy widths of the Bloch bands as a function 
of the lattice depth. For shallow lattices the band gap A2 is 
much smaller than Ai. When moving the lattice this allows 
for Landau-Zener transitions out of the 2""* band while atoms 
in the I''' band are dragged by the lattice. At Vb = 2.7hw one 
may drive a transition from the 2"^* to the S**" Bloch band 
with energy difference hQ2,5 while a transition from the first 
Bloch band is forbidden by a band gap. 



FIG. 4. (Color online) Admixtures of different eigenstates for 
different initial states during a CROT on the left control qubit 
at a: = and the central target qubit at x — tt. The Gaussian 
beam with waist a = 7r/2 at xq = O.Gvr is linearly ramped 
to a strength of 7 = 0.224 during a time of 2h/Er. After 
waiting for 12h/Er the beam is linearly ramped for 16h/Er 
to the avoided crossing with the repulsively bound state at 
7 — 0.204. After shaking the lattice resonantly to the en- 
ergy difference between |111) and |101) the gaussian beam is 
switched off in reversed order. For a gate time of 400h/Er a 
fidelity of 99.4% is reached (see text). 

quently the atom distribution of the remaining atoms by 
fluorescence imaging with single-site resolution. The re- 
moval of excited atoms is related to an evaporative cool- 
ing of the system and corresponding strategies may be 
applied. One method is to accelerate the lattice in x 
direction. For sufficiently deep lattices the large gap be- 
tween the first and second Bloch bands inhibits Landau- 
Zener transitions such that atoms in the ground state |1) 
are dragged by the lattice. However, atoms in state |0) 
may tunnel to higher Bloch bands and eventually leave 
the lattice [57| |2H] ■ The tunneling can be enhanced by 
first transferring atoms in state |0) to a higher Bloch 
band. As shown in Fig. [5] for the exemplary lattice depth 
of Vq — 2.7huj]^ one can drive transitions from the 2"^^ to 
the high-lying 5"^ Bloch band by means of shaking the 
lattice with frequency ^2,5 while transitions from the first 
Bloch band arc inhibited due to a band gap. In the 5th 
Bloch band the atoms are quasi-free and leave the lat- 
tice even during slow accelerations. Lately, a closely re- 
lated scheme of a state-dependent removal of atoms from 
an optical lattice has been used experimentally for cool- 



In conclusion, a scheme for quantum computation 
with Bosons in Id optical lattices was presented. The 
qubits are encoded in the spacial atomic wave function 
which suppresses decoherence due to fluctuating mag- 
netic fields. It was shown that by shaking the lattice one 
may drive single and controlled qubit rotations and how 
qubit dephasing can be prevented by using refocussing 
pulses. The qubit readout can be performed by removing 
atoms in excited states from the lattice and determining 
the atom distribution by fluorescence imaging. For gate 
times on the order of milliseconds fidelities above 99% 
can be reached. We believe that gate times and fidelities 
can be further improved, e.g., by using optimal control. 
It would be interesting to extend the approach to a 2d 
OL in order to reduce the number of needed swap oper- 
ations. However, considering ""Li atoms in a d = 300 nm 
lattice, already within the Id approach a factorization of 
15 would be feasible within about 50 ms with about 20 
ms needed for swapping operations. 
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